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Observational constraints on supermassive dark stars 
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ABSTRACT 

Some of the first stars could be cooler and more massive than standard stellar models 
would suggest, due to the effects of dark matter annihilation in their cores, ft has 
recently been argued that such objects may attain masses in the 10 4 -10 7 M Q range 
and that such supermassive dark stars should be within reach of the upcoming James 
Webb Space Telescope. Notwithstanding theoretical difficulties with this proposal, we 
argue here that some of these objects should also be readily detectable with both the 
Hubble Space Telescope and ground-based 8-10 m class telescopes. Existing survey 
data already place strong constraints on 10 7 M© dark stars at z w 10. We show that 
such objects must be exceedingly rare or short-lived to have avoided detection. 
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1 INTRODUCTION 



are one 



WIMPs (Weakly-Interacting Massive Particles) 
of the most promising candidates for da r k ma tter (DM; 
see e.g. Ijungman etahl Il996l ; iBergstroml I2000T ). WIMPs 
may congregate and annihilate in stellar cores, affecting 
the appearance and evolution of their host stars (e.g 
Salati fc Silklll989l; iMoskalenko fe Waill2007l; ISpolvar et al 



200a iFairbairn et al.l 
2009). 



2008; [Iocco et a l. 2008; Scot t et al 



WIMP annihilation in the dark matter halos that 
hosted galaxy and primordial star formation may have af- 
fected the ability of halo baryons to form stars there, by 
heating the gas enough to compete with radiative cooling 
l|Ascas ibar 2007|). Of particular interest is the idea that the 
first stars (Population III) might have been substantially 
affected, not only because this heating competes with the 
cooling necessary to form a star, but because it can lead to 
a new phase of stell ar evolution domina ted by annihilation 
rather than fusion (Sp olvar et al.l 2008). It h as been sug- 
geste d that this phase could be long-lived (e.g. iFreese et alj 
120081 ). but evidence in this direction is far from conclu- 
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sive l|Ripamonti et al.1 12010| ; ISivertsson fc Gondold bOloD . 
If they are long-lived , such stars could be detectable 
l|Zackrisson et alj [201oh with the upcoming James Webb 
Space Telescope (JWST), scheduled for launch in 2014. In 
this scenario, the dominant effect is gravitational contrac- 
tion of the host halo by baryonic infall a s sociated with the 
formation of the s tar (jSpolvar et al.ll2008l ; llocco et al.ll2008l ; 
IFreese et al]|2009D . Following the initial contraction, colli- 
sions of through-going WIMPs with stellar nuclei, and their 
subsequent scattering to lower-en ergy orbits, might brin 
yet more dark ma tter into the star (|Gouldlll987l ; lloccdl20oT 
IFreese et al]|200Sl ). 

Because such dark matter-powered stars (also com- 
monly dubbed "dark stars") are cooler than normal Pop- 
ulation III stars, they are able to accrete more mass than 
their canonical cousins before radiative feedback develops, 
leading to stellar masses of up to ~ 10 3 A'/ Q . To a first ap- 
proximation, the limiting stellar mass is set only by the point 
at which the star can no longer be sustained by dark matter 
annihilation; eventually, either the dark matter is depleted 
by annihilation, or the star becomes too massive for the 
annihilation to support it. In realistic models, one might 
also expect aspherical accretion, stellar winds and three- 
dimensional gas flows in the upper atmosphere to cause 
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some degree of fragmentation, probably leading to some- 
wha t smaller fina l mass es . 

iFreese et al] (|2010h have recently suggested the exis- 
tence of supermassive dark stars (SMDS), with masses of 
up to 10 Mq. As is to be expected, such mega-stars would 
be easily detectable with JWST, even at redshift z ~ 10-15. 

In order to achieve such masses, the authors postulate 
an effectively boundless well of dark matter from which a 
star might draw its power. They argue that the depletion of 
the dark matter accessible to the growing dark star can be 
wholly avoided by the existence of chaotic orbits in triaxial 
halos, as these will refill the region of phase space s hown to 
be depleted very quickly bv lSivertsson fc Gondolol |2010l ). 

However, the depleted WIMPs belong to a qualitatively 
different population than those on chaotic orbits. The de- 
pleted population consists entirely of bound WIMPs, whose 
orbits lie within the star's sphere of influence. These are 
the WIMPs that are gravitationally contracted into a dark 
star, and go on to power it. WIMPs arriving on chaotic or- 
bits come from an unbound population. Members of the un- 
bound population typ ically do not contribute signi ficantly to 
the annihilation rate (jSivertsson fc Gondoldl2010l ). because 
their occupation time in the star is extremely low (i.e. they 
are moving very quickly, and do not return to pass through 
the star again). 

It is very difficult to exchange WIMPs between the two 
populations. The only way for unbound WIMPs to repop- 
ulate the depleted region of phase space is to scatter on 
stellar nuclei, and become gravitationally bound to the star. 
Without scattering, unbound WIMPs can only significantly 
contribute to the annihilation rate if very many of them pass 
through the star per unit time. In this case, their small indi- 
vidual occupation times would be offset by the sheer number 
passing through the star. The SMDS scenario hence requires 
either an extremely high rate of centre-crossing by WIMPs 
on chaotic orbits, or a large centre-crossing rate combined 
with a very large nuclear scattering cross-section. 

WIMP-nucleon scattering cross-s ections are s t rongly 
constrained by experiment (e.g. lAbbasi et" ID [2009; 
I Ahmed et al]|2009l ), requiring a < 10 cm for scattering 
on hydrogen. Without detailed simulations of triaxial WIMP 
halos, it is difficult to comment on exactly how hard it is to 
achieve a high rate of centre-crossing by WIMPs on chaotic 
orbits. In a spherical halo, with an isotropic distribution of 
WIMP velocities outside the star's sphere of influence, and 
assuming a — 10 -38 cm 2 , the resultant conversion rate of 
unbound to bound WIM Ps is ~8 orders of magnitud e too 
small to produce SMDS (ISivertsson fc GondoloTbOloT ). For 
SMDS to be viable, this additional factor of 10 s must be 
compensated for entirely by the radial bias of chaotic or- 
bits. This would seem quite improbable, though more de- 
tail ed calculations are r equired to definitively rule it out. 

IFreese et al.l (|2010r) do not explicitly calculate the rate 
of change of dark matter annihilation in their stars, sim- 
ply assuming that the dark matter population required to 
support a given stellar mass is always present. One might 
expect substantial differences in the evolutionary history 
of such stars if the dark matter density was provided self- 
consistently. For example, any downwards perturbation in 
the dark matter annihilation rate leads to a contraction of 
the dark star. Th is leaves behind a shell of dark matter 
l|locco et al.ll2008h , further reducing the total annihilation 



rate. The only way to avoid this is for the gravitationally- 
contracted dark matter to thermalise in the protostar very 
quickly, changing its distribution from effectively flat with 
radius to strongly peaked at the centre of the star. Again, 
this requires a very larg e nuclear scattering cross-section. 

IFreese et all (feoiol ) rightly question the overall stabil- 
ity of SMDS, particularly with respect to general- relativistic 
corrections to their gravitational potentials. Regardless of 
their gravitational (in)stability though, SMDS should also 
suffer from radiative-hydrodynamic instabilities. As ex- 
pected in radiatively-supported objects, their luminosities 
all lie on the Eddington limit. This means that they prob- 
ably exhibit significant continuum-driven stellar winds, not 
to mention line-driven ones from H and He. The effect of 
including these processes in detail would almost certainly 
be to limit their m asses to something lower than found by 
IFreese et alj l|2010l ). Any upward perturbation in the dark 
matter annihilation rate would substantially accelerate this 
effect. The growth of SMDS hence requires a very delicately 
tuned, smoothly varying annihilatio n rate, in ord e r for the 
accretion to continue as proposed by IFreese et alj l|2010l ). 

As well as theoretical problems with the SMDS forma- 
tion scenario, there are also observational constraints. Here, 
we use stellar atmosphere models to estimate the apparent 
magnitudes of these dark stars. We argue that there is no 
need to wait for the JWST to observationally test the SMDS 
hypothesis, since the most massive (~ 10 7 M ) SMDS arc 
sufficiently bright to be readily detected with the Hubble 
Space Telescope (HST) or 8-10 m class telescopes on the 
ground. In fact, such SMDS would need to be exceedingly 
rare and/or short-lived to have evaded detection until now. 



2 THE APPARENT MAGNITUDES OF 
SUPERMASSIVE DARK STARS 

We use the TLU STY stellar atmosphere code 
l|Hubenv fc Land 1 19951 ) to generate spectra i n the 0.015- 
300 11 m wavelength range for the SMDS of IFreese et al.l 
(2010). To derive broadband fluxes, these model spectra 
are redshifted to z = 0-15 and convolved with the relevant 
filter transmission profiles and detector sensitivities. The 
resulting broadband fluxes are converted into apparent 
magnitudes using the luminosity distance, assuming 
H = 72 km s" 1 Mpc -1 , Q M = 0.27 and fi A = 0.73. The 
calibration is based on AB magnitude system, defined so 
that an object with a constant flux per unit frequency 
interval of 3631 Jy has zero AB magnitudes m in all filters. 
To simulate the Gunn-Peterson trough due to the opaque 
intergalactic medium during the reionization epoch, all 
fluxes are moreover set to zero at restframe w avelengths 
short ward of Lya (A < 0.1216 /im) at 2 > 6 (|Fan et al.l 
120061 ). 

The resulting broadband fluxes are valid at the sur- 
faces of these objects, and neglec ts the fact that many of 
the SMDS of IFreese et alj (|2010h also are sufficiently hot 
to photoionize the surrounding gas. The resulting Hfl re- 
gions will add emission lines and a nebular continuum to 
the observed spectra of these SMDS, thereby substantially 
boosting their fluxes in the rest-frame UV and optical (e.g. 
ISchaererll2002l) . Because of this, the fluxes we derive should 
be regarded as conservative. 



Constraints on supermassive dark stars 3 




Figure 1. The apparent - H16O AB mag nitudes of 10 5 -10 7 Mq 
SMDS from iFreese et all j2010l) as a function of redshift. The 
coloured lines represent dark star models with 10 7 Mq and T e ff = 
27000 K (thick solid blue), 10 7 Mq and T cff = 51000 K (thick 
dashed blue), 10 6 Mq and T cS = 19000 K (thick solid green), 

10 6 Mq and T eff = 51000 K (thick dashed green), 10 5 Mq and 
T cff = 14000 K (thick solid red) and 10 5 Mq and T cff = 16500 
K (thick dashed red). The thin solid red and thin dashed red 
lines (which nearly overlap in this diagram) represent the 10 5 Mq 
dark stars with T cff = 50000 K and T cS = 51000 K, respectively. 
The dashed horizontal line indicate s the detection thres hold of 
the WFC3 HUDF observations (e.g. iBouwens et al.|[2009ll . Both 

10 7 Mq models (thick blue lines) are readily detectable above the 
threshold at z = 10, and so is one of the 10 6 M models (thick 
solid green line) . The vertical lines indicate how the models would 
be shifted by a gravitational magnification of /i = 10 (short arrow) 
and fi = 100 (long arrow) from a foreground galaxy cluster. A 
lensing boost of fi = 10 would lift the remaining 10 6 Mq model 
(thick dashed green line) and two of the 10 5 Mq models (thick 
solid and dashed red lines) into the detectable range. A boost of 
fi = 100 would make all the plotted models detectable. 



In Fig. [TJ we plot the redshift evolution of the apparent 
magnitudes for 10 5 -10 7 Mq SMDS in the iJST F160W filter 
(hereafter H 160 ). At z = 10, both 10 7 Mq SMDS (thick blue 
solid and dashed lines) and one of the 10 6 Mq dark stars 
(thick solid green line) from IFreese et al l l|201oi ) are readily 
detectable above the Hiso — 28.8 mag threshold of current 
Wide Field Camera 3 (WFC3) images taken of the Hubble 
Ultra Deep Field (HUDF) in August 2009 (HST GO 11563: 
PI Illingworth) . Similar observations of fields gravitationally 
lensed by foreground lensing clusters could in principle lift 
the remaining SMDS in this diagram above the detection 
threshold. This is illustrated by the two arrows, which in- 
dicate how the lines would be shifted by magnifications of 
(j, = 10 (short arrow) and fi = 100 (long arrow). By compar- 
ison, a galaxy of similar mass would be considerably fainter 
than these SMDS. Using the lZackrisson et alJ l|200ll . 120081 ) 
population synthesis model, we predict that a galaxy with 
stellar mass 10 7 Mq, metallicity Z = 0.001 and Salpeter 
IMF (mass range 0.08-120 Mq) should have H 160 > 29 mag 
at z = 10, i.e. at least 2.5 mag fainter than a 10 7 Mq SMDS. 

To be able to separate z = 10 objects from low-redshift 
interlopers, data in more than one filter is required. Since 
the Gunn-Peterson trough renders the flux in the F110W fil- 




Figure 2. The apparent magnitudes of 10 7 Mq SMDS as a func- 
tion of redshift. Thin lines represent the T e g = 27000 K a nd th e 
thick lines the T off = 51000 K SMDS from IFreese et al.l feulCl) . 
The blue lines represent fluxes in the Jno filter and red lines 
fluxes in the Hxw filter. The vertical dashed lines mark the detec- 
tion l imits of the deepest survey data compiled by l lBouwens et al.l 
2008). The rapidly dropping J- band fluxes at z > 6 are due to 
the Gunn-Peterson trough, making such objects appear as J-band 
dropouts (Jno — Higo > 1.3) at z > 9. An similar decline is also 
seen in the Hieo band at z > 10.5. 



ter (hereafter Jno) extremely low at z « 10, objects at this 
redshift will remain undetected in Jno (and at all shorter 
wavelengths) while still being potentially detectable in Hiqo. 
Because of this effect, candidate z 10 objects can be se- 
lected from large multiband datasets by searching for J- 
band dropouts. This is demonstrated in Fig. [2] where we 
plot the red s hift e volution of the two 10 7 Mq SMDS from 
IFreese et all (|2010h in both Jno and Hi6o bands. The Jno 
band fluxes (blue lines) drop rapidly at z > 6 as the Gunn- 
Peterson trough enters this filter, making the z ~ 10 fluxes 
extremely low (implying very red J — H colours). Eventu- 
ally, the Gunn-Peterson trough also starts to affect the i/i6o 
filter, but not until z > 10.5. 

Since the most massive SMDS should be sufficiently 
bright to already be readily detectable, one might ask what 
existing deep near-IR survey data have to say about the ex- 
istence of such objects. For this purpose, we use t he the Jno 
and Hieo data compiled bv lBouwens et all l|200Sl ) from HST 
and groundbased images of the regions around the HUDF 
and the Hubble Deep Field-North. While the more recent 
WFC3 data reach slightly deeper than these, they have not 
yet been as thoroughly analysed, and the number of bona 
fide J-band d ropouts in the WFC3 images still rem ains con- 
troversia l (cf.lYan et all 120091; IBouwens et~aH 120091 ) . 

TheH 

ouwens et al.l ((2008) compilation contains a num- 
ber of survey areas with slightly different detection thresh- 
olds in the Jno and -ffi60 passbands, reaching a maximum 
depth of Jno = 28.4 and Hieo = 28.2 at 5er (dashed blue 
and red horizontal lines in Fig. [2}. By adopting a colour 
criteria of Jno — -Hieo > 1.3 for J-band dropouts, and cross- 
correlati ng potential detectio ns with data at other wave- 
lengths, IBouwens et al.l (|2008l ) reported a null detection of 
credible J-band dropouts in these images. Still, Fig. [2] in- 
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Figure 3. The formation rate of 1-2 X 10 s Mq halos per comoving 
Mpc 3 and year, as a function of redshift. The raw simulation data 
is represented by the thin line, whereas the thick line traces a 
second-degree polynomial fitted to the data. The halo formation 
rates used in the analysis are based on the latter. 



dicates that the two 10 7 Mq SMDS should appear as J- 
band dropouts according to these colours criteria, at least 
throughout the redshift range z « 9.5-10.5. This conflict 
can only be resolved if 1O 7 M SMDS are so rare or short- 
lived that no such objects are expected within the surveyed 
regions of the sky. In the following sections, we attempt to 
convert this criterion into a quantitative constraint on the 
formation rate and properties of 1O 7 M0 SMDS. 



UPPER LIMITS ON SUPERMASSIVE DARK 
STARS 



According to lFreese etaD l|20ld ). 10 7 M Q SMDS form inside 
~ 10 s Mq dark matter halos. This scenario is controversial, 
since a transition from molecular (mainly H2 and HD) to 
HI cooling is expected in this halo mass regime, supposedly 
leading to the formation of entire stellar populations ('first 
galaxies') rather than single population III stars within a 
halo. In the standard picture of chemical enrichment in the 
high-redshift Universe, short-lived population III stars in 
~ 10 5 -10 6 Mq minihalos at z > 15 moreover pollute the 
intergalactic medium with metals to such a degree that only 
a tin y fraction (< 10~ 3 in the models of IStiavelli fc Trentil 
I2009T I of the > 10 s Mq halos that form at z = 10-15 can 
be che mically pristine and able to form metal-free SMDS 
the way IFreese et all |201(t ) envision this process. The stan- 
dard scenario for chemical enrichment can be questioned, 
however, if a substantial fraction of these first population 
III in minihalos go through a long-lived dark star phase, 
since this would delay the production of pair-instability su- 
pernovae ljloccoll2009l ) and slow down cosmic chemical evo- 
lution. For the sake of argument, we put these theoretical 
problems to the side, and instead attempt to observation- 
ally constrain the fraction /smds of ~ 10 s Mq CDM halos 
that form 1O 7 M0 SMDS. If the parameter /smds is suffi- 
ciently low at z m 10-15, this would explain the lack of such 
objects in existing survey data. One should keep in mind, 
however, that both cosmic chemical evolution, the chang- 



ing intergalactic radiation field and the time-dependent halo 
merger rate may allow for a /smds which evolves within this 

redshift in terval. 

While IFreese et alJ l|201Ch consider SMDS lifetimes r of 
up to r w 1 Gyr, this parameter remains poorly constrained 
at the current time, and one may also consider values that 
are several orders of magnitudes smaller. If so, SMDS would 
only be detectable during a brief period after their forma- 
tion, hence contributing to the lack of detections in existing 
data. 

In Fig. [3l we plot the formation rates of 1-2 x 10 8 Mq 
dark matter halos as a functio n of redshift, bas ed on high- 
resolution N-body simulations ijlliev et al.ll201ol ) of the for- 
mation of high-redshift structures. We use the CubeP 3 M N- 
body c od jj] which is based on t he particle-mesh code PM- 
FAST ijMerz. Pen fc Trad l2005h and simulate a cubic co- 
moving volume of size 6.3 hT 1 Mpc using 1728 3 particles of 
mass 5.19 x 10 3 M o . To identify halos we use a spherical 
overdensity halo finder with overdensity parameter fixed to 
178 and a minimum number of particles equal to 20, i.e. the 
minimum halo mass is 1.04 x 10 5 Mq. The background cos- 
mology is based on the WMAP 5-year data combined with 
constraints from the baryonic acoustic oscillations and high- 
redshift supernovae (Qm = 0.27,0a = 0.73, = 0.7, Qb = 
0.044, cr g = 0.8, n = 0.96; iKomatsu et al J 120091 ) . 

While IFreese et al ] (|201fl ) adopt a formation redshift of 
z = 15 for their 10 7 Mq SMDS in 10 s Mq halos, it seems 
reasonable that such objects should be able to form also 
in similar halos at slightly lower or higher redshifts, albeit 
possibly with different formation probabilities /smds • In the 
following, we present constraints for two different scenarios: 
A) in which 10 7 Mq SMDS can form in halos down at least 
z ~ 10 as well, and B) in which the /smds evolves so strongly 
with redshift that 10 7 Mq SMDS effectively form only at 
2 « 15. 

At z « 10, the formation rate of 1-2 x 10 8 Mq halos is 
dn/dt « 5x 10 -9 per comoving Mpc 3 and year (Fig.[3}. This 
converts into « 580 halos formed per arcmin 2 in the redshift 
interval z = 9.5-10.5 where 10 7 Mq SMDS would be read - 
ily detectable as J-band dropouts bv lBouwens et al l |2008r i. 
This implies that, if the formation of long-lived 10 7 Mq 
SMDS within 10 s Mq halos were a common phenomenon at 
these redshifts, a survey area like the HUDF (11 arcmin 2 ) 
would light up like a Christmas tree from the glow of thou- 
sands of bright J-band dropouts. Since this is clearly not the 
case, these objects must be exceedingly rare or very short- 
lived. 

The null detections of J-band dropouts in a given survey 
area can be converted into constraints on /smds and r using 
the expression: 



/smds ^ 



At 
N6 2 t'" 



(1) 



where N is the number of 1-2 x 10 8 Mq halos forming 
per unit redshift and arcmin 2 , B 2 is the angular survey 
area in arcmin 2 and At is the cosmic age interval per unit 
redshift. At z = 10, N » 580 and At « 6.6 x 10 7 yr. 
A detailed comparison of the -H160 fluxes of the 10 7 Mq 



http://www.cita. utoronto.ca/mcdiawiki^index.php/CubePM, 
for a description of the code see also llliev etafl j200Sf) . 
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Figure 4. Upper limits on the fraction /smds of 1—2 X 10 Mq 
dark matter halos that form T e fj = 27000 K (solid line) and T e ff = 
51000 K (dashed line) 10 7 Mq dark stars at z fa 10, as a function 
of their lifetimes r. 



SMDS with the detection thresholds of the surveyed areas 
reveals that 2 = 18.5 arcmin 2 have been imaged to suffi- 
cient depth to detect the T cff = 27000 K, 10 7 Mq SMDS 
as J-band dropouts, whereas the corresponding area for the 
T cft = 51000 K, 10 7 M e SMDS is 6 2 = 3.3 arcmin 2 . From 
this we derive the two sets of upper limits on /smds as a 
function of r included in Fig. [3] 

The resulting constraints are strong: e.g. log 10 /smds =Ss 
-3.2 (-2.5) if t ~ 10 7 yr and log 10 /smds < -2.2 (-1.5) 
if t ~ 10 6 yr for the T cff = 27000 (51000) K SMDS. These 
upper limits formally apply only to the value of /smds at- 
tained at z = 10±0.5. Lifetimes in excess of At = 6.6 x 10 7 yr 
would make the constraints even stronger, but would require 
additional assumptions concerning the redshift evolution of 
/smds- 

In scenario B, where /smds is assumed to be effectively 
zero at z = 10, the existing data can still be used to set upper 
li mits on /smds at z =15 (the formation redshift assumed 
bv lFreese et~alll20ld ), provided that the SMDS forming at 
2 = 15 have sufficiently long lifetimes to survive until z = 10. 
In the adopted cosmology, this requires r > 2.1 x 10 s yr. By 
adopting a formation rate of dn/dt ~ 1 x 10~ 9 per comoving 
Mpc 3 and year for 1-2 x 10 8 Mq halos (see Fig. [3}, we arrive 
at log 10 /smds < -2.9 (-2.2) for th e T eB = 27000 (51000) 
K, 10 7 M Q SMDS from iFreese eTall (|2010l ). 



4 DISCUSSION 



As demonstrated bv lFreese etafl l|20ld ). JWST can detect 
~ 10 s Mq SMDS out to z « 10 and ~ 10 7 Mq SMDS 
out to z as 15. However, the constraints already placed on 
~ 10 7 Mq SMDS by existing data imply that the prospects 
of detecting objects in this mass range with the JWST may 
be rather bleak. At z = 15 we predict N « 30 halos with 
mass 1-2 x 10 8 M© forming per arcmin 2 and unit redshift. 
With At w 2.6 x 10 7 yr and a field of view covering 9 = 4.84 
arcmin 2 , a single JWST detection of a ~ 10 7 Mq SMDS at 



z = 15 would suggest log 10 /smds ~ —1.8 if r = 10 7 yr. 
However, this combination of /smds and r has already been 
ruled out at z = 10 (Fig. [3}. Hence, if /smds and r are 
approximately the same at z = 15 and z = 10 (as in our 
scenario A), our constraints predict that no 10 7 Mq SMDS 
will be detectable within a single JWST field at z — 15. Of 
course, JWST observations would still be highly relevant for 
dark stars at lower masses. As demonstrated in Fig.[T] HST 
observations through lensing clusters may also be able to set 
constraints on SMDS down to masses of ~ 10 s Mq. 
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